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We investigate the roto-breathers recently ob- 
served in experiments on Josephson ladders sub- 
jected to a uniform transverse bias current. We de- 
scribe the switching mechanism in which the num- 
ber of rotating junctions increases. In the region 
close to switching we find that frequency locking, pe- 
riod doubling, quasi-periodic behaviour and symme- 
try breaking all occur. This suggests that a chaotic 
dynamic occurs in the switching process. Close to 
switching the induced flux increases sharply and 
clearly plays an important role in the switching 
mechanism. We also find three critical frequencies 
which are independent of the dissipation constant 
a, provided that a is not too large. 



The recent discovery || of discrete breathers in 
Josephson ladder arrays (Fig. |l|) driven by a trans- 
verse dc bias current has shown not only that these 
localised excitations exist but that they exhibit re- 
markable behaviour, which is the subject of this 
paper. 

The Josephson ladder has been studied theoret- 
ically for many years. In the absence of a driv- 
ing current the interaction between vortices has 
been found to be exponential and this leads to 
the vortex density exhibiting a devil's staircase 
as the magnetic field is increased EL Quan- 
tum fluctuationsjs] , meta-stable states^ and in- 
ductance effects 0] have been studied as has the 
interaction of vortices with a transverse dc bias 
current^, ^. 

In the case of a Josephson ladder, a roto-breather 
is a stable group of vertical junctions rotating to- 
gether {9j — 9j w cat, see Fig. Until recently, 
the interest has focussed on the case of a sinusoidal 



bias current, where it has been shown that roto- 
breather solutions exist and that a discrete breather 
may even include a chaotic trajectory [lO|. 

Roto-breathers should also be stable in the 
experimentally easier case of a ladder array sub- 
jected to a uniform dc transverse driving current. 
The case of rotation occurring at a single vertical 
junction has been studied nume rically|l|,|l3). Two 
experimental groups |]l| ^ have now independently 
confirmed the existence of such solutions and dis- 
covered further unpredicted behaviour. In an an- 
nular ladder[|l] (Fig. |l|a) and in a linear ladder[^ 
(Fig. |l|b) breathers were observed with various 
numbers of "vertical" (i.e. radial) junctions rotat- 
ing. However the most interesting features of the 
observations were the switching between breathers 
and their symmetry breaking behaviour. Once a 
breather has been initialised it is stable, but if the 
bias current is slowly decreased then, at some crit- 
ical value of the current, the number of rotating 
junctions switches to a larger number and a new 
breather is formed. Furthermore, the new breather 
may have a different centre of symmetry from the 
old breather despite the symmetry of the ladder it- 
self and the symmetry of the driving currents about 
a single vertical junction; this is particularly re- 
markable in the case of the annular ladder where 
translational invariance should be exact. The pro- 
duction of breathers with an even number of ro- 
tating junctions is, by itself, a demonstration of 
symmetry breaking. The symmetry breaking effect 
was not mentioned by the experimenters, presum- 
ably because it could be due to experimental im- 
perfections. However we show that it should also 
arise in a perfect experiment, the apparent mech- 
anism being the occurrence of chaotic dynamics in 
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the switching region. 

First we construct a new model, similar to 
but significantly different from those proposed 
previously|l|, |l|, |lO). The current through a 
junction is determined by the RCSJ model 
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the following pair of coupled differential equations 
for each plaquette: 
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where is the change in superconducting order 
parameter 9 across the junction and A is the vec- 
tor potential. The "vertical" (i.e. radial) junc- 
tions may differ in area from the "horizontal" junc- 
tions by an anisotropy parameter rj = Ich/Icv = 
Ch/Cy = Rv/Rh where I^h (lev), Ch (C^,) and Rn 
(Rv) are, respectively, the critical current, capaci- 
tance and resistance of a horizontal (vertical) junc- 
tion. From Fig. ^ we see that there are three un- 
knowns per plaquette: 9j, 9'j and fj, where fj is the 
total flux threading the jth plaquette. To solve for 
these unknowns we construct three equations per 
plaquette as follows. The first two equations are 
obtained from current conservation at the top (in- 
ner) and bottom (outer) rails. The third equation 
is obtained by making the approximation that the 
induced flux fj — fa (where fa is the applied flux) is 
produced solely by the currents flowing around the 
immediate perimeter of the jth plaquette (Fig. ^) : 
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where [3l is an inductance parameter. This last 
equation makes the plausible assumption that the 
plaquettes are more or less square and that it is only 
the currents flowing around the plaquette that pro- 
duce the induced field. This differs from, and for 
square plaquettes is more accurate than, the com- 
mon assumption that the induced field is 
proportional to the loop (or "mesh" ) current circu- 
lating the plaquette. 

We impose the initial conditions, mimicing the 
experiments, that at i = 0, 9j = 9'^ = 0, fj = and 
Ij = for all j. This means that 9j + 9'j =0 for 
all time, i.e. there is a symmetry between the inner 
and outer rails. Using Landau gauge we then have 
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We now focus on determining whether or not our 
model exhibits the interesting switching and sym- 
metry breaking behaviour observed in the data of 
Binder et aZQ. All parameter values are chosen 
to mimic the experimental setup i.e. 77 = 0.44, 
(3l = 2.7, a = 0.07 and fa = 0. Let 
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where Ib is called the bias current. Again follow- 
ing the experiment (Fig. |^), we slowly increase 
/a while keeping Ib — until rotation starts at 
site j = 0. /a is then slowly decreased while at the 
same time increasing Ib to keep Iq constant. When 
Ib has reached the desired value it is then held 
fixed while I a is slowly reduced to zero. Finally Ib 
is slowly reduced while keeping /a = 0. Note that 
the dynamical equations, initial conditions and in- 
jected currents are exactly symmetrical about site 
j — 0. One would expect only solutions which are 
also symmetrical about j = 0. However, like the ex- 
periments, the simulations also produce breathers 
which are not symmetrical about j ~ 0. Fig. |^ 
shows how the number Nj^ of rotating junctions 
changes as the bias current Ib is slowly reduced 
from various starting values. In agreement with 
experiment, breathers with even Nji are commonly 
produced (these cannot be symmetrical about j — 
0), and Nfi switches to larger and larger values un- 
til eventually all junctions are rotating. While we 
have found that similar behaviour occurs in a pre- 
viously published model[|l^ |T^, our model gives 
considerably better agreement with the experimen- 
tally observed switching currents, thus indicating 
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the importance of inductance effects in the switch- 
ing mechanism. We beheve the origin of the sym- 
metry breaking is the occurrence of chaotic dynam- 
ics in the switching region. 

The chaotic nature of the switching is further 
suggested by the fact that when the same computer 
code which produced Fig. |^ is performed on a differ- 
ent computer (different floating point processor) we 
see significant changes in the switching behaviour 
of nearly ah trajectories although the overall pat- 
tern remains identical. The minute differences in 
the handling of floating point numbers are ampli- 
fied in the chaotic regime to produce significantly 
altered trajectories. 

Fig. H shows the voltage-current characteristics 
obtained from the same simulations used in produc- 
ing Fig. ||. At high frequencies {V = {dd~/dt) > 
4.3) most of the breathers show purely resistive be- 
haviour (i.e. V — {dd^ /dt) = IbR), the value of 
the resistance R increasing with Nj^ according to 
an expression deduced from experiment 0| : 

aR^ 1/(1 + 7^/Nr) (7) 

This is the relationship expected if the current 
through each rotating vertical junction is aV and 
the four horizontal junctions surrounding the ro- 
tating region each carry current ^rjaV (to satisfy 
Kirchoff's rule). For = 1 this expression was 
derived in ref. [|l3). We find that no such resistive 
behaviour occurs for V < 4.3. 

A typical example of a resistive breather (far 
away from any switching region) is shown in Fig. |^. 
Although the breather is not symmetric about j = 
it shows exact symmetry about the midpoint of 
the rotating region and also appears to be exactly 
periodic, the period being two revolutions of a ver- 
tical junction (we call this "period 2"). The non- 
rotating junctions oscillate, the amplitude decreas- 
ing exponentially with distance from the rotating 
region. Far away from switching the magnitude of 
the flux fj is everywhere small (< 0.1) and is lim- 
ited to the edges of the rotating region. 

Fig. ^ shows that there are at least two critical 
frequencies, {d9^ /dt) = 6.5 and {d9~ /dt) = 5.0, at 
which breathers become unstable and switching oc- 
curs. As the current is reduced the frequency falls 
until it reaches the critical value, at which point 
the breather becomes unstable and switches to a 
larger breather with a larger resistance (Eq. (R)) 



and therefore a higher frequency. The process 
then repeats as the current continues to be ramped 
down. We identify the larger of these two critical 
frequencies with that observed experimentally Q. 
Although, at fixed Is , the frequency of rotation de- 
pends on the dissipation constant a, we find that 
these two critical frequencies are more or less in- 
dependent of a, as is the critical frequency below 
which breathers cease to show constant resistance. 
This of course breaks down when the bias current 
required to achieve the upper critical frequency ex- 
ceeds unity, i.e. when a > 0.15. 

Note also that in the above simulations: (i) no 
breathers are seen for Ib > 0.85, (ii) only the sin- 
gle site breather is seen for 0.68 < Ib < 0.85 and 
(iii) all rotating junctions normally stop together 
when the current Ib is reduced below 0.17, al- 
though sometimes a new breather containing only 
a few rotating junctions is produced which is then 
destroyed when Ib falls below 0.14. 

As well as resistive period 2 breathers we have 
found frequency-locked breathers. In fact the usual 
Nji = 1 resistive breather makes a transition to 
a frequency-locked breather as the bias current is 
reduced below Ib = 0.699. At this point the 
period doubles, and the breather becomes asym- 
metric and jumps to a higher rotation frequency 
{d0-/dt) = 6.32. At Ib = 0.6895 the period dou- 
bles again and at Ib — 0.6892 it becomes quasi- 
periodic (Fig. Switching to multi-site breathers 
occurs at Ib ~ 0.6891. Such frequency locked 
breathers have not been reported experimentally 
but have probably been overlooked as they are only 
stable in narrow current intervals. 

While the maximum flux f™-"-^ is normally small 
(^jrnax ^ 0.1), it sharply increase at switching to 
jrnax ^ \ \i would appear that an important part 
of the switching mechanism is the instability caused 
by having a large flux concentrated in a small area. 
The importance of the flux in the switching process 
is further confirmed by the fact that when Ib has 
finally been reduced all the way to zero we find that 
the ladder may contain one or more stable vortex- 
antivortex pairs. 

The occurrence of frequency locking, period 
doubling, quasi-periodic behaviour and symmetry 
breaking as switching is approached strongly sug- 
gests that the origin of the symmetry breaking is 
the occurrence of chaotic dynamics in the switch- 
ing region. From the theory of chaotic dynamics 
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it is known that two or more coupled Josephson 
junctions (or, equivalently, two or more coupled 
pendula) may exhibit three types of motion: os- 
cillatory, rotational and chaotic. Chaotic motion 
arises for particular initial conditions. An anal- 
ogous situation arises in non-linear coupled lat- 
tices which also display both breathers and chaotic 
dynamics [p^ . In a large array of coupled Joseph- 
son junctions (or, analagously, an array of coupled 
pendula) we must also expect these three types of 
motion. While roto-breathers are the most charac- 
teristic stable solutions it appears that the switch- 
ing between breathers occurs only in the chaotic 
regime. Of course, the chaotic dynamics has an 
influence on the selection of breather type. Any 
small difference in the initial conditions leads to a 
completely different breather. Thus, in this chaotic 
regime, any small experimental imperfections or 
small inaccuracies in a computer simulation may 
lead to different chaotic trajectories and hence dif- 
ferent roto-breathers. Symmetry breaking arises if 
the perturbation itself breaks symmetry. 

We conclude that our model exhibits most of the 
main features (l6) of the roto-breathers recently ob- 
served in Josephson ladder arrays, sheds light on 
the switching mechanism and predicts further ob- 
servable behaviour. The occurrence of frequency 
locking, period doubling, quasi-periodic behaviour 
and symmetry breaking in the switching region sug- 
gests that switching occurs in the chaotic regime. 
We find that the maximum flux increases sharply 
in the switching region and that the flux clearly 
plays an important role in the switching mecha- 
nism. We also find two critical switching frequen- 
cies and a critical frequency below which no resis- 
tive behaviour is observed. All three frequencies 
are approximately independent of the dissipation 
constant a (for a < 0.15), 
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Figure 2: Result of smoothly ramping down the 
bias current. Each small circle marks the start of a 
new ramp. 



Figure 1: (a) An annular ladder subjected to a uni- 
form transverse bias current Ib- Each circle repre- 
sents a superconducting island. Each link between 
islands represents a Josephson junction. 9 is the 
phase of the superconducting order parameter and 
/ is the flux threading a plaquette. (b) A linear 
ladder, (c) Explanation of the notation used in 
Eq. (I). 
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Figure 3: Resistance aR (where R^V/Ib) plotted 
against voltage V — {d9~ /dt) for the data shown in 
Fig. H Note the critical switching frequencies and 
frequency locked states. 
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Figure 4: An example of a resistive breather far 
from switching. dOj /dt is plotted against t for all j. 
Note that six junctions (—2 < j < 3) are rotating 
and that the breather is not symmetrical about the 
injection site j = 0. The breather is periodic and 
exactly symmetrical about its centre. The peak 
flux (not shown) is small 0.06). 




Figure 5: The quasi-periodic breather (approxi- 
mately period 8) which occurs near Ib — 0.6892, 
just above the bias current at which switching to 
multi-site breathers occurs. It is frequency-locked 
at {d6-/dt) = 6.32. Note that d9~ /dt is nearly 
anti-symmetric about j = 0. The maximum flux 
(not shown) is large (0.5) and increases still further 
at switching. 
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